Integer Quantum Hall Transition and Random SU{N) Rotation 
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We reduce the problem of integer quantum Hall transition to a random rotation of an A''- 
dimensional vector by an su{N) algebra, where only TV specially selected generators of the algebra 
are nonzero. The group-theoretical structure revealed in this way allows us to obtain a new series 
of conservation laws for the equation describing the electron density evolution in the lowest Landau 
level. The resulting formalism is particularly well suited to numerical simulations, allowing us to 
obtain the critical exponent v numerically in a very simple way. We also suggest that if the number 
of nonzero generators is much less than N, the same model, in a certain intermediate time interval, 
describes percolating properties of a random incompressible steady two-dimensional flow. In other 
words, quantum Hall transition in a very smooth random potential inherits certain properties of 
percolation. 
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PACS numbers: 73.40.Hm, 64.60.Ak 

1. The quantum Hall transition is a delocalization 
transition of a particle moving on a two dimensional 
plane with a random potential and a strong magnetic 
field perpendicular to the plane M. Recently, an ele- 
gant approach to this transition was suggested by Sinova, 
Meden, and Girvin m. The idea of the method is to con- 
sider the quantum mechanical states of a particle belong- 
ing to the lowest Landau level, and to project the density 
operators p{x, t) onto these functions. The spatial corre- 
lation properties of the particle can then be described 
by the correlation function of these projected density 
matrices, G{x,t) = Ti{p{x,t)p{0,0)) or by its average, 
{G{x,t)). In this formula, "Tr" corresponds to summing 
over all the states of the lowest Landau level, and the 
angular brackets denote averaging over the random po- 
tential. 

In a certain sense, G{x, t) is the probability of tran- 
sition of a quantum particle from the origin to point x 
in the course of time t. This statement would be ex- 
actly true if we allowed the particle to travel over all 
the quantum states, i.e., not only those belonging to the 
lowest Landau level. Not having the full system of wave 
functions does not allow us to localize the particle at a 
distance smaller than the magnetic length l'^ — hc/{eB). 
For example, the initial condition for the G function have 
the form G{x,0) = Aex-p{—x'^/2P), with some normal- 
ization constant A. However, for strong enough magnetic 
fields, the magnetic length I is mush smaller than the size 
of the system and particles can be adequately localized. 

The projected density operators obey the Schrodinger 
equation 






[H,p], 



(1) 



whose Hamiltonian is just the random potential V{x) 
projected onto the lowest Landau level, H = 
Jd X V{x)p{x). The correlation function G satisfies 
an analogous equation, which due to remarkably sim- 
ple commutation relations between the density opera- 



tors p, can be written in the Fourier space {k is a two- 
dimensional vector, k = {ki,k2)) 0: 

ih—G{k, t)= I d^q 2i sin (-kxq\ V{k - q) x 

G{q,t). (2) 



exp 



--{k^-k-q) 



For a simple, "from the first principles," derivation of this 
equation we refer the reader to f^. All the information 
about the quantum Hall transition (QHT) is contained 
in this equation, which we use as a starting point in the 
present work. Unfortunately, it is not obvious how to 
solve this equation. 

In this letter, we are going to demonstrate the hid- 
den group-theoretical structure of Eq. (g) . We show that 
it admits a series of conserved integrals of motion that 
have not been known before. We reduce the problem 
of QHT to a finite-dimensional problem of a random ro- 
tation by banded su{N) matrices with the band width 
~ VjV. We also show that when the width of the band is 
much smaller than \/N, such a model describes classical 
percolation. 

Let us start with the classical limit 1^0. In this case, 
Eq. (0) is simplified considerably. 



d_ 
di 



G{x,t) =eijdiV{x)djG{x,t). 



(3) 



A formal solution to (H) can be written at once [pj. It is 



G{x,t) = 6{x — x{t)) where 



dxj 

dt 



eijdjV{x), the equa- 



tions which were thoroughly studied in M. From here it 
is easy to deduce that we are trying to describe a par- 
ticle which percolates along the equipotential lines of a 
random potential. Obviously that does not capture the 
physics of quantum Hall transition. 

The reason why a solution to (0) was so easy to write 
down is contained in its infinitely many integrals of mo- 
tion. Indeed, it is not hard to check that any expression 



of the form 



= / d'x G'^V" 



(4) 



is conserved along the solutions of (^ . The conservation 
of Iio and In is simply a consequence of the conserva- 
tion of probability and energy. However, higher order 
integrals have a less obvious meaning. Now, in a quite 
remarkable way, these integrals of motion do not get de- 
stroyed as one goes back to the original equations (||) , at 
least for integer m and n. They only get slightly modi- 
fied. 

To see that, let us put the particle on a torus. The 
wave functions of the lowest Landau level can be written 
explicitly in the Landau gauge. 



V'a(a;,y) 



y expl 2tt{x + iy){N'm + a) 
{Nm + af 



N 



-ttNx^ 



(5) 



Here < x, j/ < 1, iV denotes the number of flux quanta 
through the torus and a goes from to iV — 1, label- 
ing the N states on the torus in the lowest Landau 
level. These wave functions describe the electron local- 
ized along a narrow strip around the line x = a/N. No- 
tice that the magnetic length I is now automatically cho- 
sen in the form 



I' 



1 



2ttN 



(6) 



It is now a matter of simple calculations to project the 
density operators P(ki,k2) — exp (27ri (fcia; -I- fc2j/)) onto 
the lowest Landau level on the torus. In these notations, 
fci and fc2 are integer numbers. As a result, the density 
operators are now N x N matrices in the basis of states 
(h) which can be written in the following form. 

Consider a pair of unitary unimodular N x N matrices: 



/ = diag(l,e, . 



N-l 



(7) 



(8) 



where /i is a cyclic permutation matrix, and e = 
exp(27rz/iV). These matrices have the following prop- 
erties: hf = efh, h^ — f^ ~ 1. Now introduce the 
matrices 



^(fei,M=e''''^'/''^'' 



(9) 



where {ki,k2) 7^ (0,0). The matrices ^(fcj^fe^) are peri- 
odic in fci and fc2 with period N up to coefficient ±1. 
They can be chosen as a basis for the su{N) algebra. 



The product and the commutator of two such operators 
have the following form. 



TTl 



LqLp = exp ( — 9 X p ) Lq+p, 



[Lq, Lp] = 2i sin ( —qxp] Lq+p. 



(10) 
(11) 



One can easily write down the explicit expression for the 
matrix elements of these matrices: 



fr \ _ ^fcl'c2/2+fei(Q-l) r I 



modw 



(12) 



Then the density operator can be expressed in terms of 
these matrices as 



PiMM) = exp 



' ^2N ^ ^) -^(fci.fca)- (1^) 



To be specific, (|13| ) gives the density operator in the 
Schrodinger representation, as opposed to (n^) where the 
Heisenberg representation was assumed. From this point 
on, we will understand p as a time-independent density 
operator. 

It is not difficult to check that the density operator 
written in this form indeed satisfies the commutation re- 
lations discussed, for example, in H with the magnetic 
length chosen according to (g). 

The projected Hamiltonian takes the form 



H ^ 



Y^ V{-ki,-k2)p{kuk2)- 



(14) 



feifeo 



Since the Fourier components of the potential V{x) are 
random, we see that the Hamiltonian is none other than 
a random su{N) matrix. This completes the projection 
to the lowest Landau level on the torus. At this point we 
can easily rederive (0) simply by writing p as a matrix 
and commuting it with the Hamiltonian H as prescribed 
in (U) using the explicit matrix definitions (|l^). Finally, 
let us introduce the matrix G — ^j. j, G{k,t)p(kx,k2)- 
Here as a consequence of choosing ki and fc2 on the torus 
as integer numbers, we denote k = 27r(fci, /C2). As a con- 
sequence of (0), it satisfies the equation identical to (1^), 



ih^G= [H,G]. 
dt ^ ^ 



(15) 



Through all these manipulations we succeeded in re- 
ducing the original problem of a particle on an infinite 
plane in a random potential to a problem with a finite 
number of degrees of freedom. That allows us immedi- 
ately to write down the generalizations of the integrals 
of motion (Eh . Indeed, since a trace of a commutator of 
finite matrices is equal to zero. 



/„„ = Tr (g"F" 



(16) 



is conserved. These integrals can be expressed in terms 
of G(fc, t) and V{k) using explicit expressions for the ma- 
trices p. 

The first few of the integrals in (p6|) coincide with their 
"classical" counterparts in (|j). For example, /lo and In 
are equal to /lo and In and are still the probability and 
energy conservation, respectively. Higher order integrals 
of motion become increasingly more complicated. For 
example, 



ho = Yl dik) g{s) g{r-k- s) 



X exp 



N 



{k X s + k X r + s X r) + TTiNrir2 



(17) 



where k — 27r(fci,fc2), s = 2Tr(si, S2), r = 2TrN(ri,r2), 
g{k) = G'(fc)exp[-7r(fc^ + fc|) /(2iV)], and the summa- 
tion is going over integer fci, fc2, si, S2, ri, and r2. 

Despite being rather complicated, this expression re- 
duces to its classical counterpart, /30 — J cPx G^{x), in 
the limit I -^ {N -^ cxj). It is of course possible to 
show, after some algebra, that it is indeed conserved un- 
der the time evolution (^. A simple mathematical rea- 
son lies behind this: in the limit A^ — > 00, the operators p 
become the generators of the group of volume preserving 
diffeomorphisms on a torus MM ■ Such a group represents 
motion of an incompressible fluid, which is precisely the 
meaning of (0). 

2. We are going to demonstrate now that all the fea- 
tures of QHT are preserved in this picture. Consider a 
particle placed in one of the states M). That means, it 
is localized along the x direction and is extended along 
the y direction. Wait some time t and measure the aver- 
age square of the displacement of the particle along the 
X direction. 
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i-'m^Tn 



N 
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(eR 



Ht 



a/3 



{a-Pf 



(18) 



where ()q./3 denotes the matrix element of the matrix in- 
side the brackets. We expect that, according to (|J^, at 
large enough times but before the finite size of the system 
is reached, the deviation ( x'^{t) ) behaves as 



(2:2)^^1-1/(2-). 



(19) 



The critical exponent u is to be determined, and is be- 
lieved to be close to 7/3 . 

The Hamiltonian (nj) appears to be a linear combina- 
tion of su{N) matrices with random coefficients. How- 
ever, if it were indeed a random su{N) matrix with prob- 
ability distribution invariant under SU{N) rotations, it 
would lead to a particle instantaneously hopping all over 
the torus. In fact, such a Hamiltonian would be unphysi- 
cal. As a consequence of the strong exponential suppres- 
sion of high Fourier modes in the density operators (13), 



the SU{N) rotation invariant Hamiltonian corresponds 
to the random potential V{x) varying considerably at 
distances much smaller than the magnetic length. In 
fact, in the continuum limit it would lead to the poten- 
tial having Fourier harmonics growing exponentially at 
large fc, as follows from the analysis of ([l^). 

A more physical setting involves the potential whose 
Fourier harmonics at least do not grow at large k. While 
it is possible to choose a particular random V{x) satisfy- 
ing this property and compute the Hamiltonian dl4), it 
is not necessary to do so. Instead, we can just point 
out that the exponential in (|l3|) strongly suppresses 
high Fourier harmonics of the random potential, such 
that fci,fc2 > n = ^IN j-K. Therefore, we can simply 
choose the Hamiltonian to be 



A;i,fc2=0 



^(fcl,/C2)) 



(20) 



where ^(/ci^fcj) are random independent variables with 
equal mean square values ( w^ ) . The crucial part of (|20| ) 
is that the Hamiltonian is a linear combination of only A^ 
generators of the su{N) algebra, out of total A^^ gen- 
erators, with random coefficients. From the explicit 
form ( [12|) of matrices L(j~-^i^^^ we see that the random ma- 
trix H has nonzero elements only on the diagonal strip 
of the width ~ N^^^. Within the strip, only the ma- 
trix elements no farther than ^ N^^^ from each other 
along a given diagonal are correlated with each other. 
Such banded random matrix theory, which reproduces the 
physics of the QHT, was considered in a similar context 
in , although the Hamiltonian ( pO| ) was not introduced 
there. 

This framework provides a very convenient setting for 
numerical calculation of v. Below we present the nu- 
merical simulations for matrix 1000x1000. The Hamil- 
tonian has been chosen in the form (pO[). The ran- 
dom coefficients V(^ki.k2) have thus been generated in the 
square j/cij, |A:2| < n = y^2N/TT, and the real and imagi- 
nary parts of each Fourier mode ^(fcjj;^) — ''^i-k -k ] ^^^^ 
chosen randomly, independently, and uniformly from the 
interval [—0.5, 0.5]. Instead of calculating the matrix ex- 
ponent in (nSh directly, we have simulated the equation 



Za — iHapZp, 



(21) 



where Za is a wave function in the representation of the 
states on the torus. It is easy to check that in this rep- 
resentation, 



g{ki,k2) = Zl [L(^ki.k2)\al3^i 



/3> 



(22) 



where g{k) is introduced in (O). 

After the random potential H had been generated, and 
the initial distribution had been chosen in the form 



Z,, 



1, if a = OfQ, 
0, otherwise, 



(23) 



the dispersion {x^{t) ) = Y^^ia — aoY\Za\^ was calcu- 
lated as a function of time. The same calculation was 
performed for all the initial positions of the particle, 
ao = Ij • ■ • I -^i a-iid the average was taken over all such 
realizations. The result is shown in Fig. (Q). As we see, 
we closely reproduce the universally accepted value of v 
without much difficulty. 




FIG. 1. Log((a;^)) is plotted as a function of log(t) for 
iV^lOOO, n = 15. The light line has the slope 0.79 which 
corresponds to v = 2.38. 



3. As was already discussed, a naive I -^ limit 
leads to the percolation picture of QHT. However, I -^ 
is the same as TV ^ oo. It should indeed be possi- 
ble to reproduce the percolation behavior if the scale of 
the random potential is much larger than the magnetic 
length Z, but much smaller than the size of the system. 
In other words, we need to keep only such modes in (llj) 
for which 1 <C n ^ N^/^. The percolating behavior 
holds until the width of the diagonal strip of the ma- 
trix exp{iHt) becomes equal to N^^^. After that, the 
regime changes to QHT. The square of this width is 
given by (ttq), and therefore grows as t^^^^^^'^ \ with the 
percolation v* equal to 4/3. The crossover time (mix- 
ing time) is thus equal, in physical units, to Votm/fi ~ 
(iV/n2)2'^*/(2i'*-i) _ (;^/;)4^V(2>^*-i)^ where lo is the 

scale of the random potential and Vq is its typical ampli- 
tude. The percolation regime is valid for t < tm- 

Our simulations of Eq. (p^) show that for n < 
(2A^/7r)^/^, the behavior of x^ coincides with ( [l9| ) at t -^ 
oo, but follows some intermediate asymptotics before 
that. We expect these asymptotics to be identical to 



those describing steady percolating flow of an incom- 
pressible two-dimensional fluid, i.e., v* = 4/3. This 
seems to agree with our numerics, but more extensive 
simulations are required f^. In the other limit, n 3> 
y/2N/TT, the behavior becomes diffusive, (a^) ~ t. 

In real experiments, the random potential length is of- 
ten much larger than the magnetic length I. It would be 
interesting to devise an experiment which would probe 
the intermediate percolation asymptotics, perhaps by 
looking at finite frequency conductivity. 

In conclusion, we have presented a model that re- 
duces the problem of quantum Hall transition to a finite- 
dimensional problem of random rotation by su{N) ma- 
trices. As a consequence, the equation (pi) describing 
quantum Hall transition, admits a series of integrals of 
motion. The random rotation matrix H is quite arbi- 
trary, except for the dependence on the parameter N 
through the exponential cutoff in (|l3|). This rather cru- 
cial dependence amounts to vanishing of all the matrix 
elements but those belonging to the diagonal strip of the 
width n ^ N^/^. Within this framework, the quantum 
Hall transition exponent v can easily be evaluated numer- 
ically. Attempts to change n independently of N lead to 
other models, e.g. those describing diffusion and percola- 
tion in a steady incompressible 2D flow. That leads to a 
prediction that in very smooth potentials quantum Hall 
transition should exhibit some properties of percolation. 
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